Phase synchronization between collective rhythms of globally coupled oscillator groups: 

noiseless non-identical case 

Yoji Kawainur4j 
Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology, Yokohama 236-0001, Japan 

Hiroya Nakao 

Department of Physics, Kyoto University, Kyoto 606-8502, Japan and 

JST, CREST, Kyoto 606-8502, Japan 

Kensuke Aral 
Brain Science Institute, RIKEN, Wako 351-0198, Japan 

Hiroshi Kori 

("^ ' Division of Advanced Sciences, Ochadai Academic Production, 

^Nj , Ochanomizu University, Tokyo 112-8610, Japan and 

PRESTO, Japan Science and Technology Agency, Kawaguchi 332-0012, Japan 

Yoshiki Kuramoto 

0^ ' Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan and 

Cn , Institute for Integrated Cell-Material Sciences, Kyoto University, Kyoto 606-8501, Japan 

(Dated: July 26, 2010) 

Cj ' Phase synchronization between collective oscillations exhibited by two weakly interacting groups 

,^rfH ' of non-identical phase oscillators with internal and external global sinusoidal coupling of the groups 

^J , is analyzed theoretically. Coupled amplitude equations describing the collective oscillations of the 

f^ ■ oscillator groups are obtained by using the Ott-Antonsen ansatz, and then coupled phase equations 

for the collective oscillations are derived by phase reduction of the amplitude equations. The col- 
lective phase coupling function, which determines the dynamics of macroscopic phase differences 
between the groups, is calculated analytically. It is demonstrated that the groups can exhibit 
effective anti-phase collective synchronization even if the microscopic external coupling between in- 
dividual oscillator pairs belonging to different groups is in-phase, and similarly effective in-phase 
collective synchronization in spite of microscopic anti-phase external coupling between the groups. 
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O . I. INTRODUCTION 



Assemblies of coupled dynamical elements are abundant in Nature. Macroscopic collective oscillations typically 
emerge in such systems through mutual synchronization of the individual elements [lH6|. Based on coupled phase 
oscillator models, theoretical investigations of the origin and the nature of collective oscillations have been carried out. 
C^ ' Globally coupled phase oscillators, a representative class of such models, have been particularly well analyzed |7H14|. 
The appearance of an experimental system of coupled electrochemical oscillators [15l - [l9| has worked so powerfully 
in accelerating the study of collective dynamics of globally coupled oscillators. Recently, Ott and Antonsen [20, l2l[ 
proposed a remarkable mathematical ansatz for the analytical treatment of coupled phase oscillators in the continuum 
limit, which is applicable to models with global sinusoidal coupling. Since then, various applications |2a - [3l| and 
extensions [33, [33] (see also Refs. [3J,[3a]) of the Ott-Antonsen ansatz have been rapidly developed. 

When two or more groups of dynamical elements exhibiting collective oscillations interact with each other, synchro- 
nization among those collective oscillations may naturally be expected. In Refs. |36l44C) |. two interacting groups of 
globally coupled phase oscillators have been studied and mutual entrainment between the groups have been reported. 
Description of the system in those works was essentially based on microscopic phases of the individual oscillators, and 
macroscopic properties of the collective synchronization were also investigated on the microscopic footing. However, 
it should be more convenient and beneficial if one can describe the collective oscillations at the macroscopic level 
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using appropriate macroscopic variables in a closed form. 

In this paper, using the collective phase of each oscillator group as the macrovariable |4l| - |44| . we formulate a 
theory of synchronization between two interacting groups of globally coupled oscillators closed at the macroscopic 
level, based on the Ott-Antonsen ansatz that gives a low-dimensional description of phase oscillators with global 
sinusoidal coupling, as well as on the standard phase reduction theory for limit-cycle oscillators. We analytically 
derive the collective phase coupling function, which determines the macroscopic dynamics of the interacting groups, 
and illustrate several representative cases of collective phase synchronization between the groups. In particular, 
we reveal a counter-intuitive phenomenon in which the macroscopic collective phase difference between two groups 
becomes anti-phase in spite of microscopic in-phase external coupling between individual pairs of oscillators between 
the groups, and also the opposite phenomenon, namely, in-phase synchronization between oscillator groups with 
microscopic anti-phase external coupling. 

In Ref. [43| , we considered a similar problem, namely, collective phase synchronization between two groups of globally 
coupled phase oscillators. The crucial difference is that we treat deterministic noiseless non-identical phase oscillators 
in the present paper, whereas we analyzed stochastic noisy identical phase oscillators in Ref. [43|. Though these two 
cases may look similar, they are physically different systems and require distinct mathematical approaches; here we 
rely on the Ott-Antonsen ansatz, whereas we used center-manifold reduction as well as phase reduction to nonlinear 
Fokker-Planck equations in Ref. [43]. As we will show, we still find very similar transitions between effective in-phase 
and anti-phase synchronization in both cases. This implies that the two systems have similar effective low-dimensional 
dynamics, despite the fact that the two models are originally defined in completely different high-dimensional phase 
spaces. 

The organization of this paper is as follows. In Sec. |IT1 we introduce a model of weakly interacting groups of 
globally coupled phase oscillators and illustrate both effective anti-phase and in-phase synchronization of collective 
oscillations by numerical simulations. In Sec. IIIIl we develop a theory based on the macroscopic phase description of 
collective oscillations that clarifies whether the phase coupling between collective oscillations is effectively in-phase or 
anti-phase. In Sec. IIV[ we illustrate several representative cases of the collective phase coupling function obtained in 
Sec, mil and reexamine the numerical results in Sec.|TTl Concluding remarks will be given in the final section. 

II. COLLECTIVE PHASE SYNCHRONIZATION OF OSCILLATOR GROUPS 

We consider two interacting groups of globally coupled non-identical phase oscillators described by the following 
model: 
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for j = I,--- ,A^ and (c, r) = (1,2) or (2,1), where (pf {t) is the phase of the j-th oscillator in the cr-th group 
consisting of N oscillators. The second term on the right-hand side represents internal coupling between oscillators 
within the same group, and the last term gives external coupling between oscillators that belong to different groups. 
This phase model can be derived from two interacting groups of globally coupled limit-cycle oscillators near the Hopf 
bifurcation point using the phase reduction method 2] under weak coupling conditions. Phase oscillators with global 
sinusoidal coupling as given in Eq. ([T]) can be experimentally realized in electrochemical oscillator systems [15l-tl9]. 

The internal coupling is specified by the parameters K > Q and |a| < 7r/2, where K determines the coupling 
intensity and a gives the coupling phase shift. These parameter values correspond to in-phase (attractive) internal 
coupling. Similarly, the external coupling is specified by the parameters J > and |/3| < tt. The external coupling 
can be either in-phase (attractive) (|/3| < 7r/2) or anti-phase (repulsive) (|/3| > 7r/2) 0. The characteristic magnitude 
of the weak external coupling is given by a small parameter e > 0. 

The natural frequency ojj is assumed to be drawn from the Lorentzian distribution with central value wo and 
dispersion 7, 



We define a parameter 



sH^Jt ^r—-,- (2) 

^=^r^— , (3) 

K cos a 

which is the ratio of the frequency dispersion 7 to the attracting component K cos a of the internal coupling function. 
When the external coupling is absent, i.e., e = 0, Eq. ([T]) describes two independent oscillator groups [a — 1,2), each 
of which exhibits collective oscillations under the condition < 77 < 1/2 [Sj. 



Introducing a complex order parameter A'-'^^t) with modulus i?(°')(i) and phase 0*^'^^(i) through 
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for each group, we can rewrite Eq. ^ as 

0^"^ = cj, - KR^"^ sin ((/)5") - e(") + a) - e Ji?(") sin (^J"^ - 9^") + /?) . (5) 

The order parameter R^'^^ quantifies the degree of synchronization in each group (0 < R^°''> < 1), and 0^"^^ gives the 
macroscopic collective phase of the group. 

Focusing on weakly coupled collective oscillations, we carried out numerical simulations of Eq. ([5]) with Eq. ^ 
under the following conditions: Without loss of generality, we can assume K — J ^ 1 and ujq = 0. The external 
coupling was assumed to be much weaker than the internal coupling, e = 0.01. We set the frequency dispersion 
7 — {K cos a) /A and the internal coupling phase shift a — Stt/S, so that the parameter ij was given by 77 = 1/4. The 
number of oscillators in each group was N — 4096, which was sufficiently large to observe clear collective oscillations. 

Figure [ija) shows typical evolution of the collective phase difference, |8*^^-' — O*^^-*], with in-phase condition for 
microscopic external coupling, /? = Stt/S. Two groups of phase oscillators exhibiting collective oscillations were 
separately prepared with their collective phases being almost equal, and these states were used as the initial condition. 
Despite the in-phase external coupling condition for individual pairs of oscillators in different groups, the collective 
phase difference became anti-phase (jG'^-* — G'^-*! = tt) after some time. Thus, Fig. [TJa) implies that effective anti- 
phase coupling between collective oscillations is realized. In contrast, Fig.[ljb) shows effective in-phase synchronization 
between collective oscillations (jO^^-' — 0*^^^! — 0) with microscopic anti-phase external coupling, /? — — Stt/S. 

Snapshots of the phase oscillators after the collective phase difference has reached the asymptotic value in Fig. [T] 
are displayed in Fig. [2l The oscillators are sorted in increasing order of their natural frequencies. The coherent 
segment represents phase-locked oscillators within each group and scattered points correspond to drifting oscillators. 
The coherent, phase- locked segment of each group is not centered about the middle oscillator because the internal 
coupling phase shift a is non-zero. In Fig. [2][a), the two distributions of the oscillators are shifted by tt, indicating 
anti-phase synchronization between the groups. In contrast, the two distributions almost overlap in Fig. ^b), i.e., 
they are in-phase synchronized. Note that drifting oscillators from different groups do not synchronize with each 
other, in other words, the collective phase synchronization between the groups is not due to complete synchronization 
of individual oscillators at the microscopic level. 

III. DERIVATION OF THE COLLECTIVE PHASE COUPLING FUNCTION 

We now derive collective phase equations describing the interacting oscillator groups via the amplitude equations 
obtained by using the Ott-Antonsen ansatz |20l . l2l| . We analytically determine the collective phase coupling function 
and its type, specifically, whether it is in-phase or anti-phase. 

Using the order parameters defined in Eq. Qj, we can rewrite Eq. ([T|) as 

0f ' - c., - - (A^^^ e^^f' e- - ^(-) e-^^'f' e^ ) - - (a(-) e^^f e^^ - A*-) e'^^^"' e"'^) , (6) 

where A^"' is the complex conjugate of A^*^'. In the continuum limit, N — >■ 00, we can obtain the following continuity 
equation for each a [20, [2l[ : 
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which describes dynamics of the probability density function /('^^((/<, w, t) of the phase and the frequency. Here, 
f^'^\(l),uj,t) satisfies normalization conditions 

dcu/(^)(0,c^,t) = l, (8) 






= 5 (w) , 


Jo 



and the complex order parameter A^"^' is now defined by 



^('^) (t) = 



for cr = 1,2. 

We now apply the Ott-Antonsen ansatz [20, [2l| 
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to the continuity equation ([7]), which replaces all the Fourier coefficients of f'^"\4>^ui^t) by integer powers of the 
complex variable a^'^'ioj^t). This ansatz leads to a two-dimensional representation of the infinite-dimensional partial 
differential equation when the frequency distribution g{uj) is Lorentzian. It has been shown that the above restricted 
functional form of f'''^\(j),u},t) yields asymptotically correct dynamics of the complex order parameter A^'^\t) |21| . 
By substituting this expression into Eq. ([7]) , we can derive the following equation for the complex dynamical variable 
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where (ct, r) = (1,2) or (2, 1). Moreover, in the case of the Lorentzian frequency distribution, Eq. ©, the complex 
order parameter A'^'^\t) can simply be expressed by a'^'^\uj,t) as 



A('^) it) = 
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du e^'^/^'") (0, u,t)= / du a^'^'l (w, t) g (w) = d^"^ {u = ujq - i-f, t) 



(12) 



by performing a complex contour integral in the lower-half complex plane of w, where w = wq ^ *7 gives the pole of 
the integrand with the Lorentzian g{ijj) [20|, |2l|. Therefore, we arrive at the following coupled amplitude equation for 
the complex order parameter A^'^i{t) in a closed form: 



A^"^ = {^i + in^) A^"'' - g 
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for [a, r) — (1, 2) or (2, 1), where the parameters are defined by 
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Note that Eq. (fT3|l describes two coupled Stuart-Landau oscillators, each of which (i.e., A ^ (m + i^c)A — g\A\'^A) 
represents collective oscillations of the respective oscillator group. Also, note that it is valid for the whole parameter 
region of the system, not only near the synchronization transition points of each oscillator group, r) = 1/2 {p, — Q). 
This is in sharp contrast to the conventional center-manifold reduction method ^ that leads to similar coupled 
amplitude equations, but which is valid only near the bifurcation points. 

As mentioned above, each oscillator group exhibits collective oscillations when r/ < 1/2 (/i > 0). Correspondingly, 
in the absence of the external coupling, e = 0, Eq. (J13p has a circular limit-cycle solution Aq{Q) on the complex plane, 
whose analytical expression can explicitly be given by 
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Lj^ — K sin a + 7 tan a. 



(15) 



where Rq and Q, represent the amplitude and the frequency of the collective oscillation in the model ^ with e = 0, 
respectively (see Appendix. \^. The right Floquet eigenvector of the limit-cycle solution Ao(8) associated with the 
zero eigenvalue is given by C/o(0) = dAo{&)/dQ, namely. 
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and the corresponding left zero Floquet eigenvector at each Q can be taken as 

c/*(e) 
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Taken together, they satisfy the normalization condition 
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= 1. 



(18) 



Though the above quantities are expressed by complex numbers for the sake of convenience in analytical calculations 
performed below, they are exactly the same as the known results for the Stuart-Landau oscillator [0]. 

Now let us introduce weak external coupling, i.e., we assume that e takes small positive values and treat the last 
term of the amplitude equation (J13p as a perturbation. Using the phase reduction method [2], we can obtain the 
collective phase dynamics of the amplitude equation (J13p by projecting it onto the unperturbed limit-cycle orbit as 



e^-^) = Re lux (e^'^A A^'^'^ 
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where we approximated A*-°'' by the unperturbed solution ^40(9^°^') and used Re[C/g (0)Ao(8)] = fl. We have thus 
obtained the following coupled collective phase equation from Eq. ([T^ : 



(20) 



for (ct, r) = (1,2) or (2,1), where the collective phase coupling function is given by 

r (e(") - e(")) = Re t;* (e^-') |jAo(eM) -d(Ao (e(")))'Ao(e("))| . (21) 

Similarly, we can also derive the collective phase sensitivity function [ij (see Appendix. |B|) . 

By inserting the expressions of Eqs. ([TH) . ([15]), and ([T7|) into Eq. (pij) . the collective phase coupling function F{Q) 
is obtained in the sinusoidal form, 

F{Q)^-psm{Q + S), (22) 

where the parameters p and S are respectively the modulus and the argument of a complex number given by 
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{(1 — ry) cos/3 — rytanasin/3} -I- i{(l — ry) sin/3 -I- 77 tan a cos /3} 



(23) 



This formula is the main result of the present paper. It determines the collective phase coupling function, Eq. (|22p in 
the collective phase equation (PU]) . which is derived from Eq. ^ via the complex amplitude equation (fT^. The type 
of coupling is found from the following quantity: 



pcosS = j|(l — 77) cos/3 — ?7tanasin/3}, 



(24) 



where pcosd > represents the in-phase coupling and pcosd < gives the anti-phase coupling. Reflecting the 
symmetry of Eq. ([1]), Eq. (|24p is symmetric about the origin in the Q!-/3 plane. 

IV. REPRESENTATIVE CASES OF COLLECTIVE PHASE COUPLING FUNCTIONS 

We here illustrate five representative cases of the collective phase coupling function obtained in Sec. IIIIl which 
correspond to several special sets of the parameters, i.e., the ratio rj given by Eq. (|3]), the phase shift a of the internal 
coupling function, and the phase shift /3 of the external coupling function. We then reexamine the results of our 
numerical simulation in Sec. |lTl 

(i) The first case is 77 = 0, which implies that all oscillators are identical, i.e., 7 = 0. In this case, the oscillators in 
the same group become completely phase synchronized due to the in-phase internal coupling, so that the maximum 
amplitude of collective oscillations is realized, namely, -Ro = 1- Inserting 77 = into Eq. (l23l) . we obtain the following 
result: 



77 = 0, 
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(25) 



which says that the parameters of the collective phase coupling function are identical to those of the microscopic 
external phase coupling function, so that the types of the effective coupling between the groups and the external 



coupling between individual oscillators coincide. The same result for the completely phase synchronized case has been 
obtained in different ways |42h44| . Note that the above results are independent of the value of the internal coupling 
phase shift a, so that a does not affect the collective phase coupling function at all. 

(ii) The second case is the limit 77 — >■ 1/2, which indicates that each oscillator group is exactly at the onset of 
collective oscillations, i.e., Rq — >■ 0. Inserting 77 = 1/2 into Eq. (|23p . we obtain the following result: 

^7^^, pe'^ = — ^e^("+^), (26) 

2 2 cos a 

which yields the real part pcos^ — (J/2){cos/3 — tanasin/3}. Thus, the microscopic internal coupling parameter 
a most significantly affects the parameters of the collective phase coupling function, in contrast to the case (i). 
Depending on the values of a and /3, the types of the effective coupling between the groups and the external coupling 
between individual oscillators can be either the same or opposite. 

(iii) The third case is a = 0, which gives an antisymmetric (odd) internal coupling function between individual 
oscillators. In this case, r] — j/K and Rq = -/I — 2?7. Inserting a — into Eq. (|23p. we obtain the following result: 

a = 0, pe'^ = {l-7j)Je''^. (27) 

Thus, the type of the collective phase coupling function is solely determined by the microscopic external coupling 
phase shift /3. Further, the collective and microscopic external coupling functions are of the same type. Similar 
scenarios have been encountered in different models |42h44| . 

(iv) The fourth cases correspond to special values of the microscopic external coupling phase shift /3, which give 
symmetric (even) or antisymmetric (odd) external coupling functions. Inserting (3 = (in-phase), ±7r (anti-phase), 
and ±7r/2 (marginal) into Eq. (P5|) . we obtain the following results: 

/3 = 0, pe^'' = J[(l-?7)+i?7tana], (28) 

/3 = ±7r, pe^^ = J[-{l-i])-iT^teiiia], (29) 
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/3 = ±-, pe'" = J[Tr?tana±z(l-77)]. (30) 

For antisymmetric (odd) microscopic external coupling functions, i.e., for /? = and ±7r, the type of the collective 
phase coupling is not affected by the internal coupling phase shift a, because a does not appear in the real part pcosS. 
In contrast, for the symmetric (even) microscopic external coupling, i.e., /3 = ±7r/2, the type of the collective phase 
coupling function is solely determined by the internal coupling parameter a. The types of effective coupling between 
groups and the external coupling between individual oscillators coincide when /3 = 0, ±7r and also when /3 = ±7r/2 
and a = 0. 

(v) The fifth case is /3 = a, namely, the case that the external and the internal couplings have the same coupling 
phase shift. Since we assume \a\ < 7r/2, the value of /3 should also be in this range. Inserting P = a into Eq. ([23]), we 
obtain the following result: 



/3 = a, pe^ = J 
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(31) 



In this case, pcosi5 depends on the internal coupling phase shift a and on the ratio rj. Therefore, effective anti-phase 
collective coupling can be realized when cos^ a < 77, in spite of the microscopic in-phase external coupling (|a| < 7r/2). 
When 77 — ?► 1/2, the effective anti-phase collective coupling is realized for cos^ a < 1/2, i.e., 7r/4 < |a| < 7r/2. 

Now, let us reexamine the case with ry = 1/4, which we used in the numerical simulations displayed in Fig. [T] The 
types of the collective phase coupling function of Eq. (|22|) is shown in Fig. [3] on the a-/3 parameter plane, where 
the solid curves represent boundaries between the in-phase and anti-phase regimes. The curves show the marginal 
condition pcosS = 0, which is determined from Eq. (|24|) . Two sets of parameter values used in generating Fig. [1] 
are plotted in Fig. [31 As can be seen, the set of parameters corresponding to Fig. [TJa) is in the effective anti-phase 
regime, whereas that corresponding to Fig. [Ijb) is in the effective in-phase regime. Therefore, the theory developed 
in Sec. [Til] successfully explains the numerical results displayed in Fig. [T] 

V. CONCLUDING REMARKS 

Three cases of macroscopic phase descriptions for collective oscillations exhibited by coupled phase oscillator systems 
have been established for (i) phase coherent states in nonlocally coupled noisy identical oscillators, (ii) fully phase- 
locked states in networks of coupled noiseless non-identical oscillators, and (iii) partially phase-locked states in globally 



coupled noiseless non-identical oscillators. Here, the case (ii) can be fully analyzed from the viewpoint of dynamical 
systems, while the other cases (i) and (iii) necessitate statistical treatments. 

The collective phase dynamics of the case (i) was established in Refs. [4ll - l43l |. where the collective phase equation 
was derived for the first time. In this case, it is essential to derive a nonlinear Fokker-Planck equation from coupled 
Langevin phase equations by using the mean-field theory [2j], which is applicable for nonlocal coupling as well as for 
global coupling in a large population of identical oscillators with independent noise. Applying the phase reduction 
method to the nonlinear Fokker-Planck equation, we can derive the collective phase equation. Furthermore, using the 
center-manifold reduction method in addition to the phase reduction method, a detailed analysis can be performed 
near the onset of collective oscillations via the supercritical Hopf bifurcation. 

The collective phase description for the case (ii) was formulated in Ref . |4J| . In this case, we can systematically treat 
any system size, connectivity, heterogeneity in the coupling, and nonuniform external forcing, as long as the oscillators 
exhibit fully phase-locked collective oscillations. In particular, the Jacobi matrix of the collectively oscillating solution 
takes the form of the Laplacian matrix encountered in graph theory [4^, |4a] , so that several analytical results can be 
obtained by using the matrix tree theorem [43, |4g . There exist several studies related to this case [4^, Ho] (see also 
Rcfs. [5T,'52l). 

The present paper provides a tractable example of the collective phase description for the case (iii) . The keystone 
in our analysis is the Ott-Antonsen ansatz [20, l21| , which is unfortunately limited to the case with global sinusoidal 
coupling and Lorentzian frequency distributions, but which yields analytically tractable coupled Stuart-Landau equa- 
tions for the complex order parameters. By virtue of the circular symmetry of the limit-cycle, we could explicitly 
calculate the collective phase coupling function between the groups. 

However, a general framework for collective phase reduction in the case (iii) is still missing. It would be necessary to 
derive a continuity equation, such as the nonlinear Fokker-Planck equation, similarly to the case (i), which can easily 
be done. However, the fundamental difficulty in applying the phase reduction method to the continuity equation in 
this case lies in the fact that the zero eigenvalue corresponding to the collective phase mode may not be isolated, but 
immersed in the continuous spectrum on the imaginary axis, as implied by the linear stability analysis [53|- In the 
formulations of the cases (i) and (ii), it is critically important that the zero eigenvalue corresponding to the collective 
phase mode is isolated. In the present study, reduction of the infinite-dimensional phase space to a finite-dimensional 
manifold by using the Ott-Antonsen ansatz yielded an isolated zero eigenvalue corresponding to the collective phase 
mode. But it is an open problem at this point how to extend the present analysis to more general phase coupling 
functions with higher harmonic terms and to more general frequency distributions of the oscillators. 

In conclusion, we have established an analytically tractable example of the collective phase description of globally 
coupled non-identical phase oscillators. We have found that the type of the collective phase coupling function can be 
different from that of microscopic external coupling function, and clarified the relation between them by systematically 
deriving the collective phase equation from the microscopic phase equations. The collective phase reduction would 
serve as a powerful method in analyzing meta-groups of coupled oscillators comprised of multiple interacting groups. 

Appendix A: Self-consistent theory of collective oscillations 

As a validation of our arguments based on the Ott-Antonsen ansatz [20, UM i '^® compare the limit-cycle solution 
given in Eq. ([THj) with the result obtained by a conventional self-consistent theory [2|. As is well known, the self- 
consistent equation [2, [3, |3] for the order parameter amplitude R of Eq. ^ with e = 0, is given in the following form 
(see also Refs. [54l. ISSj): 
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where contributions from both coherent and incoherent parts are expressed in a single formula |56|-[58j. There is a 
unique eigenvalue il of the collective frequency for which the self-consistent equation (JA1|) of the order parameter 
amplitude R admits a solution. For the Lorentzian distribution of Eq. ^, we can analytically solve the self-consistent 
equation (jAip as follows. No such calculation seems to have been carried out so far. 
The Lorentzian distribution of Eq. ([2]) can be expressed by 

.9 (^) = - Yt V^4^ ^ n ■ (^2) 

TT [[uj - ujq) + ij\ [{lo - Wo) - lj\ 

Taking the upper half-plane as the contour of integration for Eq. (jAip , we obtain 

Re-^VT^^ + ^^, ,^ i-^-^)^^\ (A3) 



This equation can be transformed into 

i?2e*" - 2iRz - e"'" = 0, (A4) 

which is equivalent to the following simultaneous equations: 
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R^ sin a ^^^ '- + sin a = 0. (A6) 
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Solving these equations, the amplitude and frequency of the collective oscillation can be respectively obtained as 

i?2 = l--^-, (A7) 

K cos a 

D, ^ uji^ — K sin a + 7 tan a, (AS) 

which coincide with the results given in Eq. (J15p . 

Appendix B: Derivation of the collective phase sensitivity function 

We consider a group of globally coupled non-identical phase oscillators subject to common weak external forcing, 
ep(t), described by the following equation: 

K ^ 
h = ^J - ]^ H sin {<t^3 - <l>k + a)+eZ {(l)j)p{t), (Bl) 

fe=i 

where the microscopic phase sensitivity function [2] is taken as 

Z(0) =-sin(/). (B2) 

The Ott-Antonsen ansatz is applicable also in this global sinusoidal case [20, l21[. Therefore, we can derive the 
amplitude equation for the complex order parameter A(t) in the following form: 

A = {n + iilc) A-g \A\^ A + ei (1 - A^) p{t). (B3) 

By applying the phase reduction method, the collective phase equation is obtained as 

e = r! + eC(e)p(<), (B4) 

where we assumed that the external forcing is sufficiently weak. The collective phase sensitivity function is given by 



C (6) = Re 



UU^)\{l-{A,{Q)f} 



which is sinusoidal in form. In general, the collective phase sensitivity function C(0) differs from the microscopic 
phase sensitivity function Z[(j)). When the oscillators become completely phase synchronized, 77 = 0, i.e., Rq — \, the 
collective phase sensitivity function coincides with the microscopic one, C,{Q) = Z{Q). Near the onset of collective 
osciUations, t] < 1/2, i.e., /^ > 0, the amplitude of the collective phase sensitivity function increases as C(0) = 
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FIG. 1: (Color online) Time evolution of collective phase difference |B'^' — 0'^'|. (a) Effective anti-phase collective synchroniza- 
tion with microscopic in-phase external coupling, /3 = Stt/S. (b) Effective in-phase collective synchronization with microscopic 
anti-phase external coupling, /3 — — 57r/8. The other parameters are K — J = 1, e — 0.01, ujq — 0, 'y — (K cos a) /A, and 
a = 3tt/8. The number of oscillators in each group is A'' = 4096. 
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FIG. 2: (Color online) Snapshots of the asymptotic states of the oscillators in Fig. [T] The oscillators in each group are sorted 
in increasing order of their natural frequencies. Only one in every two oscillators is plotted. Open circles (o) and plus signs (-I-) 
indicate oscillator in group (1) and in group (2), respectively. The collective frequency of each group is il = (—3/4) sin(37r/8), 
which differs from the central frequency wq ~ 0. (a) Effective anti-phase coupling with microscopic in-phase coupling, /3 = Stt/S. 
(b) Effective in-phase coupling with microscopic anti-phase coupling, /3 — — Svr/S. 
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FIG. 3: (Color online) A diagram showing whether the collective phase coupling function between the groups is in-phase or 
anti-phase on the a-/3 plane for a G (— 7r/2,7r/2), /3 £ [— 7r,7r], and t] = 1/4. The solid curves are theoretically determined from 
Eq. (|24|) . i.e., by pcosS = 0. The filled circle (•) indicates a = /3 = Svr/S corresponding to Figs. [Ha) and[2ja), and the cross 
(x) indicates a = Stt/S and /3 = —'o-k/?, corresponding to Figs. [2b) and[2jb). 



